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^ ■ Abstract 

We investigate the Baxter equation for the Heisenberg spin model corresponding to 
a generalized BFKL equation describing composite states of n Reggeized gluons in the 
multi-color limit of QCD. The Sklyanin approach is used to find an unitary transformation 
from the impact parameter representation to the representation in which the wave function 
factorizes as a product of Baxter functions and a pseudo-vacuum state. We show that 
the solution of the Baxter equation is a meromorphic function with poles (A — ir)~^ n ~ l > 
(r = 0, 1, ...) and that the intercept for the composite Reggeon states is expressed through 
the behavior of the Baxter function around the pole at A = i. The absence of pole 
singularities in the two dimensional A -plane for the bilinear combination of holomorphic 
and anti-holomorphic Baxter functions leads to the quantization of the integrals of motion 
because the holomorphic energy should be the same for all independent Baxter functions. 

W Work supported partly by INTAS grants 1997-31696, 2000-366, CRDF grant RP1-2108, by 
NATO and by the Russian Fund of Fundamental Investigations 

1 Introduction 

In the leading logarithmic approximation (LLA) of perturbative QCD the Reggeons (reggeized 
gluons) move in the two-dimensional impact parameter plane ~J? and interact pairwise |T], ^]. 

*Laboratoire Associe au CNRS UMR 7589. 



To unitarize the QCD scattering amplitudes at high energies one should take into account the 
multi-Reggeon exchanges in the t-channel. The composite states of the reggeized gluons satisfy 
a Schrodinger-like equation ||. 

The Reggeon Hamiltonian in the infinite color limit N c — ► oo takes a simple form and can 
be written as follows ||, 

H= l -{h + h*) , [h,h*] = 0, (1) 

where the holomorphic and anti-holomorphic Hamiltonians 

n n 

h = J2 Kk+i , h * = Yl K,k+x > ( 2 ) 

are expressed in terms of the pair BFKL operator [1], |4[]: 

h k ,k+i = logpfc + logpfc+i H (log p k ,k+i) Pk H (logp fcifc+ i)p fc+ i + 27 . (3) 

Pfc Pk+i 

Here p fc , fc +i = Pk - Pk+i , Pk = i ^ , P* k = i ^r , and 7 = --0(1) is the Euler-Mascheroni 
constant. 

In this context the Pomeron is a compound state of two reggeized gluons and the Odderon 
is constructed from three reggeized gluons. 

The operator h is invariant under the Mobius transformations with generators: 

n 

M=J2M k ; M 3 k = Pk d k , M k ~ = d k , M+ = -p 2 k d k . 
fe=i 

The Casimir operator of this group is 

n 

M 2 = -J2pld t d r . (4) 

Kr 

The Hamiltonian h describes the integrable XXX spin model with the spins being the 
generators M k of the Mobius group |J. The integrals of motion of this model are generated 
by the transfer matrix which is the trace of the monodromy matrix satisfying the Yang-Baxter 
equation |5|]. Therefore, the Quantum Inverse Scattering Method |6|, 0] can be applied to find 
an algebraic solution of the Schrodinger equation. 

The pair hamiltonian (|3]) can be obtained from the fundamental monodromy matrix asso- 
ciated to the XXX Heisenberg spin model || ||, [|. Notice that the local operators p k , p k act in 
an infinite dimensional Hilbert space whereas the spin operators in the usual Heisenberg model 
are finite dimensional matrices both for integer or half plus an integer spin. 

The auxiliary L-operator for the Heisenberg spin model with s = — 1 is given below || [!| |10| 

L k (u)=( U + P \ Pk0 ~ Pk ). (5) 

V Pk p k0 u-p k p k0 j 

where po is the coordinate of the composite state. 



The auxiliary monodromy matrix for this model can be parametrized as follows 

T(u) = L n {u)L n -i(u) ...I*{u)=( £M *M ) , (6) 

where n is the number of the reggeized gluons. The transfer matrix is the trace of the mon- 
odromy matrix 



t(u) = A(u) + D(u) = J2Q. 



n-j 



3 u 



3=0 

where 

H>in>—>ij 

The eigenvalues A(u) of t(u) take the form 

A(u) = 2u n + q 2 u n - 2 + q 3 u n - 3 + . . . + q n . 

where qj, 2 < j < n are eigenvalues of the integrals of motion Qj ||. In particular, q 2 = 
—m(m — 1) is the eigenvalue of the holomorphic Casimir operator (|J) and m is the conformal 
weight. 

The operator C{u) annihilates the pseudovacuum state for the s = — 1 model ||, 

n 

c(u) n = o , n = n p;o 2 • (8) 

r=l 

The operators B(u) can be obtained directly from eq.@. We find for n = 2 and n — 3, 
£(™= 2 )(u) = -u(p 1 +p 2 )+Pip 2 pi2 (9) 

J B (n=3) (u) = -U 2 {pi +P2+P3) +U (pi p 2 Pl2 + Pi P3 Pl3 + P2 P3 P23) - Pi P2 P3 Pl2 P23 

For arbitrary n one obtains, 

n— 1 n 

fi W M = -^&fcM n-1 ~ fc wh ere b = P = Y,Pi > 61 = - H Pi Pj Pij , 

k=0 i=\ l<*<i<« 

C*2 = / ,, Pil Pl2 P«3 P«l*2 P*2*3 > 

l<«l<i2<«3<1 



"2 V -U / v Pu P«2 • • • Pit Pii+i Pi\ii Pt2h ■ ■ ■ Piik+i i 

l<il<i2<---<ii+i<n 



b n -l = (-1)™ 1 P1P2 ••• Pn Pl2 P23 ■ ■■ Pn-l,n • (10) 

The operators B with different spectral parameters commute 

[B(u) , B(v)} = 

and therefore, one can write them in factorized form as a product of the operator zeros A^ of 
B(u): 

n— 1 



B(u) = -P JJ (u - AfcJ , AfcuAfe, 



fe=i 



At., P 







following E. K. Sklyanin JT2 



The wave function describing composite states of reggeized gluons in the holomorphic impact 
parameter space p can be written as follows [Q (see also [|]) 

n 

l/j(Pl, P2, ■ ■ ■ , Pn\ Po) = Q(Ai)Q(A 2 ) ... Q(A n _i) JIPk) 2 > PrO=Pr-pO- (H) 

r=l 



The function Q (A) satisfies the Baxter equation ||T3 



A(A) Q (A) = (A + i) n Q (A + *) + (A - «) n Q (A - i) . (12) 

For the Odderon case, the dependence of the energy from the eigenvalues of the integrals of 
motion has been found with the use of the Baxter equation [113] and of the duality symmetry 
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In this paper we systematically develop the construction of composite Reggeon states using 
the Baxter-Sklyanin (BS) representation, in which the operator zeros Xk of B(u) are diagonal. 
The matrix elements relating the momentum and BS representations obey solvable ordinary 
differential equations for n = 2, 3. These matrix elements are elementary functions for the 
pomeron case and hypergeometric functions for the odderon case. In the BS representation the 
wave function of the composite state is written as a product of the Baxter functions and the 
pseudo-vacuum state. 

For the pomeron, we provide general formulas for the Baxter function valid in the whole 
complex A plane and study its analytic properties. It turns out that the most efficient way to 
solve the Baxter equation in the present context is to use the pole expansions (Mittag-Loffier). 

We show that the Pomeron wave function has no singularities on the real axis as a function 
of a = Re A and hence it can be normalized. This corresponds to the single- valuedness condition 
in the coordinate representation. 

We derive also the analytic Bethe Ansatz equations and construct the Baxter function as 
an infinite product of Bethe Ansatz roots. 

The solution Q(u) of the Baxter equation for the general n-reggeon case is constructed as an 
infinite sum over poles of the orders from 1 up to n — 1. Their residues satisfy simple recurence 
relations. It is shown, that the quantization condition for the integrals of motion follows from 
the condition of the cancellation of the pole singularities in the two dimensional A -plane for 



the bilinear combination of holomorphic and anti-holomorphic Baxter functions Q(X) and the 
physical requirement, that all Baxter functions with the same integrals of motion yield the 
same energy. 

For the odderon, we explicitly construct the BS representation and investigate the properties 
of the odderon wave functions in this representation. The completeness and orthogonality 
relations for these functions are discussed. 

We derive new formulas for the eigenvalues of the Reggeon hamiltonian written through 
the Baxter function. These formulas generalize the result for the Pomeron to any number of 
Reggeons. The energy turns to be expressed in terms of the behavior of the Baxter function 
near its poles at A = % which are present for arbitrary n. 

The BS representation promises to be an appropriate starting point to find new composite 
Reggeon states for n > 3. In particular, it will be interesting to generalize the Odderon solution 



constructed in ref. [15] to the case of many Reggeons. 



2 BS representation for the wave function 

In order to solve the Baxter equation, one should fix the class of functions in which the solution 
is searched. The case of integer conformal weight m has been considered in refs.[[|, [10, I8f . It 
was assumed there that the solutions were entire functions with the asymptotics 

g(A)~A m -", 

but such functions do not exist for physical values of the conformal weights m. 

We want to find the conditions which should be satisfied by the solutions of the Baxter 
equation from the known information about the eigenfunctions $ of the Schrodinger equation 
in the two-dimensional impact parameter space Jf [Q, 16 . For this purpose we perform an 
unitary transformation of the wave function $ to the BS representation in which the operator 
B{u) is diagonal. 

To begin with, let us go to the momentum representation (with removed gluon propagators): 



n „ n 

iK^) 2 / n 



d 2 p k 



2vr 



exp(ipt ■ p k0 ) 



$m,m(Pl 5 P~2i . ..,P„;p )- (13) 



r=\ " fc=l 

Here $ m ~ (pi , p~2, . . . , p^; po) is the wave function of the composite state in the two-dimensional 
impact parameter space ~~j?. It belongs to the principal series of the unitary representations of 
the Mobius group and is an eigenfunction of its Casimir operators 



M 2 $ ~ = m(m — 1)$ ~ , (m*) $ ~=m(m — 1)$ 
Here 



m,m ' 



1 n _ 1 n 

m = — h iv H — , m — — Y iv 

2 2 ' 2 2 



are conformal weights (the quantities v and n are correspondingly real and integer numbers for 
the principal series of the unitary representations). The Casimir operators of the Mobius group 
are given by eq.(Q). 

For example, for the pomeron and odderon we have respectively |fj. ||, 



(\ m / % \ i 



,(-') 



where is a function of the anharmonic ratio, 

P12P30 
x == 

P10P32 ' 
Due to the identity 

n p n n—1 k , p\ n 

J2PkPk = —J2Pk + Z Pk,k+\ Z [Pr ) ' P = E Pk , 

fc=l n fe=l fc=l r=l V n/ fc=l 

we can express the quantity B(u) in momentum representation in terms of P, p& — — and the 
operators 

■( 9 d \ . d , , n 

It is convenient to introduce the new independent variables 

P =y,Pk, *i = In — , t 2 = In — , . . . , £„,_i = In . (14) 

£~J -P - Pi -P - Pi - P2 Pn. 

The quantities £& take their values in a strip of the complex plane 

—00 < Retk < 00 , — 7r < Imtk < vr . 
There is a helpful representation for the operators b\ besides that given by eq. (JiT)|). Namely, 

bi = (-lY E ( Z>J ( E P n )...( E Pn+i) flPw s+ i, (15) 

l<il<i 2 <— <h<n \ri=l / \r 2 =ii+l / \r-; +1 =i ; +l / s=l 



which is related with the duality transformation JTT] consisting of the cyclic permutation 



Pk — > Pk,k+1 — > Pk+l 



and the transposition of the operator multiplication. Note, that with the use of the duality 
symmetry the wave function in the momentum space (|l3|) for P = is proportional to the 



same function in the coordinate space |lT[. Furthermore, the function for arbitrary P can be 
obtained by an appropriate Mobius transformation. 

In the variables t%, . . . , £ n -i the matrix element B(u) takes the following form 



n-l 

B(u) = -J2hu n - 1 - k 

k=0 



where the operators bk are given by 



fe-i 



h = P £ 

l<h<l2<—<lk<n—l r=l 



Y[ (l - e* v ~Vi 



k Q 

S =l Wis 



This representation can be obtained from eq.fll5|) for hi taking into account the formulas: 



% 0_ = AE k r= lPr)(P-j: k r=1 Pr) 

l dt k P 0(£* =1 Pr 



P 



Pk,k+1 



'r-> 



1-e 



tu-t. 



PiUsllP. 



(El= + iPs)(P-Ei=iP, 



In particular, for n = 2, 3, 4 we obtain 



>(2) 



S^fu 



-P u + z 



A' 
0*1 



^(3), 



u 



-P 



u2+lu (k + i)-^- e " t2 )kk 



B^\u) = -P 
where 



r=l l<h<l 2 <3 k=l 



, (16) 



^i/ 2 — hi ~ ti 2 , d r 



_d_ 
dt. 



Since in the momentum representation the norm of the wave function is given by 

n d 2 Pr 



*. 



n 



1 \Pr 



*m,m(Pl. P2. ■■■,?£ 



we obtain after extracting the factor S 2 (P— J2k=i Pk) from \& m ~ in the new variables ti, £*; . . . ; £ n -i, C-i 



. n— 1 n— 2 



* 



m,m 



/n d2 *riil i - et *" t * +i 

J 1 1 



-2 



* 



m,m 



r=l s=l 

The operators -B(w) are symmetric B = B l with respect to this norm with the weight 

ra-2 



n 

s=l 



= t s — is+i 



-2 



The eigenvalues of the operator zeroes A& and X* k of B(u) in the holomorphic and anti- 
holomorphic space have the form, 



Afc = <?k + i 



N, 



\l = a k -i 



.N k 



where a k is real and N k is integer. The unitary transformation between t- and A- representations 
conserves the norm of the wave function 



*. 



Let us introduce the kernel U-, 



n-l ( . +O0 +oo \ 

n / *>r E 

r =l y-co N r =-ooJ 



*. 



ti, t 2 , ■ ■ ■ , t n -i) for the unitary transformation be- 



Ai,...,A n _i 

tween the t- and A- representations. It satisfies the eigenvalue equations 



n— 1 



5(u) £/■ 



Ai,...,A n . 



■ci) = -pno*-*;)^ 



fe=i 

n— 1 



Ai,...,A n _i 



(ti , . . . ,t„_i) , 



B{ u yu-> — ->(ti,..., t n _i; 

Ai,...,A n _i 
and the orthogonality relations 

•"- 1 rf 2 t fc "" 2 



-p nfa-**)* 7 



fc=i 



Ai,...,An-l 



tn— 1 1 



(17) 



5 ( 2 ^ 



J-j- | 1 _ e *r-t. + l|-2 ^ 



[/* 



7? 



r=l 



Ai,..,A„_i Al-X-i 



71 — 1 



s F n[^-o^ 



iV fc ,AT' 



t=l 



Here, A£, A& stands for the eigenvalues of the operators \ k and A£, respectively and the symbol 
Sp means the sum over all possible permutations ri, ..., r n _i of the indices 1, 2, ..., n — 1. 

To construct the kernel with the correct normalization, let us take into account, that the 
Kronecker and Dirac 5-functions appears in the right hand side of the above equation as a 
result of the integration over the region 



U - U » 1, U - U » 1 £„ i - Co >> 1 



t-3 — L 2 



, . . . , L n —\ vn—2 



corresponding to one of the two possibilities 

Pi «P2 « ••• «Pn 

or 

Pi » P2 » ■■■ » Pn- 

In this region the operator B(u) simplifies as 

n— 1 

B( U )^-P H(u + id k ) 

k=i 



(19) 



and therefore the kernel for the unitary transformation corresponds to the Fourier transforma- 
tion 



U, 



Ai,...,An-l 

where $ are some phases. 



t% , ti , . . . , t n -i) — ► 2 2 Ep e 



^V^ «**(An,-,Ar^i) 



n-1 



exp[iY,(tkK k + t*kK)} 



fe=i 



3 BS representation for the pseudo-vacuum state 

The wave function of the pseudo-vacuum state in the momentum representation is 



/n 
n 



7T \pko\ 



exp(ipt ■ p~k Q ) 



]J \p k \ 2 \n\p k \ 
fe=i 



(20) 



where we subtracted from the distribution \p\ its singular part [19 . This function is an 
eigenfunction of the transfer matrix with the eigenvalue 

A(u) = (u + i) n + (u-i) n . 

In particular, its Mobius conformal weights are 

m = rh = n 

and q k = for odd values of k. Notice that the Baxter function for the pseudo-vacuum state is 
u- independent. 

It is important to know the wave function of the pseudo- vacuum state in the BS representa- 



tion. We see from eq.(16) that Ai in the new variables for the Pomeron state takes the simple 
form 

Oti 



and the change of the basis results in 



<p 1 p 2 |P,Ai,Ai >= 



P 



P1P2 



;s* ; (i)"' -" 



Pl ~ P2) 



(21) 
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One can obtain from eqs. ( |20| ) and fl2~T|) the Pomeron pseudo-vacuum wave function in the 
new variables as 



•»(?,£) = \p\ l j<e P (^p (r^) "'hpI 2 '"I 1 



It can be written as follows 



d d 



$ u (^,A 1 ) = |Priim— — /^p|p| 



-2ia 



P 



2«/i 



y(i-p) 



p 



where Ai = a + iN/2. The integral is calculated in a closed form with the use of the anti-Wick 
rotation (see next section): 



$°(P, A 



4tt(-1) 7V |A| 2 lim . l .. . 



(22) 



The fact that the wave function for the pomeron pseudo- vacuum state turns out to be divergent 
is connected with the fact that such state having the weights m = m = 2 is outside the space of 
physical states. However it is important to normalize correctly the Baxter function. The natural 
regularization of the Heisenberg model can be provided by changing the spin representation 
s — » — 1 + e without losing its integrability. It would lead in particular to the modification of the 
pseudo-vacuum state (|20"D and to the convergence of integrals after their analytic continuation. 
The above result for <&°(~P, A ) can be obtained in a simpler way by taking into account, 
that in the integral transformation to the A-representation the large momenta p dominate 



®°(P, T 



-1^ 



\P\ 



j2 — i— -«4i 2 i |2 / 1 \N | 7->|4 I \ |2 

a pe p e p In \p\ = (— 1) \P\ |A| c, 



where the leading divergent contribution to c does not depend on A. It corresponds to the 
following simplification of the operator B(u) for p 3> P, 



B^i 



U) = - P ( u - P p t i / 



In the case of three particles we have for large pi, P2, P3 with fixed P, 



S (3), 



U 







d 2 



U +2iU di- y \dy> 



ST 

dt 2 



(23) 



-i 



t = t 1 +t 2 = P [p^ - p 3 ') < 1 , y = t 1 -t 2 = P[p- 1 
Its eigenvalues and eigenfunctions for small t and y are 

B^( u ) = -P(u-Xl)(u-X* 2 ), 
tp\ tx *(t,y) = eito+W [1 - K A* y In y) , ^(t,y) = e^+^y 

10 



Pz 1 ) « 1 • 



(24) 



Imposing the property of the single-valuedness we can write the transition amplitude in the 
two-dimensional space: 



u T 1 X2 {tu h 



c ^ _> e l[( A i + ^)*+(Ai+A 2 )t*] ( y I y 

Ai,A 2 UiA 2 AJA^ 



-In \y\ 



where the constant c^ v* ^ or ^ ne normalized function U is calculated below (see (|124j) ). 
Ai, A2 
However, the pseudo-vacuum state and other solutions with integer conformal weights do 

not belong to the space of the physical states. Therefore their unitary transformation to the 

BS representation should be special and the normalization of the unitary transformation could 

include only the integration over the large momenta. In this case it is more natural to write 

for the kernel of this transformation the following expression 



u r r (t u t 2 



Ai,A2 



where the constant c?f> — > does not depend on Ai j2 . 
Ai,A2 
Let us now consider the n-Reggeon case. Again, the large momenta pi, 1 < k < n 

presumably dominate the pseudo- vacuum wave function $°(i^, A ) when expressed as integral 

transform of eq.(^). This large momenta regime corresponds to small t^, 1 < k < n — 1 

according to eq.([U 



Retk 



P 



Z^r=l Pi 



+ 



P 



,Lr=lPr 



Imtk = 7T . 



Since the operator B^ n \u) for small Retk and Imtk — n contains more derivatives dk than 
factors t nr2 compensating them, we obtain in this regime for an arbitrary number n of reggeons 



B (n) { 



U 



-P 



n-l 



U ' 



+ u n ~ l %{n -1)% + ... 



-P 



n-l 



n— 1 



u 



.n-2 



U 



Ea- + ... 



i=i 



where t = t\ + . . . + £ n -i- 

Therefore, the transformation kernels for n reggeons have the large-p behavior similar to 
the case of n = 2 and 3, 



n-l 



U, 



AlvAn-l 

where the constant c 



>{ti , . . . , t n -\ 



d'l 



En — 1 * v~^ n — 1 

1 = 1 l + .Z^3 = l 



[] {e 1 ^^^^^^^ 



Ai,...,An-l k=l 

ps is fixed by the normalization condition with the integration 



Ai,...,A n _i 

over the region of large momenta. 

Thus, the pseudo-vacuum state in the BS representation can be written as follows 



<S>°(P,Xu---^i)~\P\ 2n £ 



Ai,...,An-l 



X 
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n-l 



ra-2 



n 



V n_1 A* y™ _1 A 

n u^^r^^H'*) ni^-wr 2 n w 2 i*k 

fe=l \ / r=l m=l 

ro — 1 , * v—*™~l 



I p\2n 



Ai,...,A n _i 



n-l 

n 

fe=i 



d 2 Pfc e 



Tl — 1 \ v-^fc v— \k I 



IT ln IP" 



m=l 



Using dimensional arguments we can write the result of the integration for the pseudo- 
vacuum wave function as 



$°(i ? „At ) ...,A^)~|P| 2n c p - 



Ai,...,An-l 



n-l 
s=l 



2(n-l) 



(25) 



up to a divergent A-independent factor which can be regularized by changing the value of the 
Heisenberg spins s — > — 1 + e. 

Thus, providing that only large momenta are essential in the unitary transformation, the 
pseudovacuum state in the BS representation is expressed in terms of the normalization constant 

2(n-l) 



Ai, 



,A n -i 



times 



Ss=i A. 



For the pomeron c 



As 



it was argued above, for the kernel U-: 



•(*! 



-r» does not depend on A [see eq.(ETT)] 
A 

,.,...., t n -i), describing the transition between 
Ai,...,A n _i 

the momentum and BS representations for the pseudo-vacuum wave function, it is natural to 

take into account in the normalization condition only the contribution from large momenta 

(presumably this is valid also for all states with integer conformal weights m and m). We 

obtain in this way, 

c£> ► = 1 (26) 

Ai,...,A n _i 

up to some A^ independent factor. 

The above arguments are in accordance with the Sklyanin theory in which the pseudo- 
vacuum state in the wave function is considered as a multiplier allowing to write the other 
multiplier as a product of the Baxter functions. Because, as it will be shown below, the 
Baxter function for the n-reggeon composite state contains the poles Q(X) ~ (A — ir) - ^ 1 ' 
for r = 0,1,..., it is natural to expect that the wave function of the pseudo- vacuum state 
cancels some of these poles. Moreover, since from each solution of the Baxter equation we can 
obtain other solutions multiplying it by factors sinh (2% A) , this symmetry should appear as a 
possibility to multiply the pseudo-vacuum state by such factors. Generally, the pseudo-vacuum 
state is not symmetric under the permutation of the parameters Ai, A2, ..., A n _i. In order to 
find the Hamiltonian in the Baxter-Sklyanin representation we show below that in the region 
where A n _i, A*_ x — > i and A s — ► (s = 1, 2, ..., n — 2) the holomorphic wave function has only a 
single pole at A n _i = i. Therefore, in order to agree with the Baxter representation we should 
substitute 



lim 

A n — 1, A n _ 



c JJN 

0,. 



,o,A 



sinh n - 2 (27rA n _ 1 ) sinh n - 2 (27r A n _ 



(27) 



n-l 
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4 BS wave function for the Pomeron 

The wave function of the Pomeron in the momentum representation is given by, 

2 'd*p k 



*m,m(Pt P*) = II (^) 2 / IT 
r=l J k=\ 



2lY 



expiipkpM,) 



Pl2 \ / Pl2 ] 

P10P20/ VPl0P20/ 



:28i 



and corresponds to the contribution of a triangle diagram. According to Appendix A the matrix 
element for this triangle diagram in the momentum space is 



where 



*m.m(Pl> Pa) = C'n 



c„ 



d 2 A; 



k(P-k) 




fc*(P* - k*) 


pt-k 




p\ - k* 


jrn-m m ~ 


r(i- 


m) 



m— 1 



(29) 



2 m+m+4 ^3 Y(m) 



Note, that the above integral over k is convergent at the singular points of the integrand 
and at the infinity providing that m and m correspond to the principal series of the unitary 
representations 

1 n _ 1 n 

m = — \- iu -\ — , m = — V iv . 

2 2 ' 2 2 



Terms proportional to 5^ (pi) and S^^pZ) []T7| can be neglected here since they are multi- 
plied by pi 2 P2 2 . 

One can analytically continue it for other values of m and m. In particular, for m = m —>■ 2, 
corresponding to the pseudo-vacuum state, we obtain that the leading contributions from two 
regions: small \pi — k\ 



2 , k*(p\+p* 2 -k* 



d l k 



p\ — k* 



, m L ( k(pi+p 2 - k) \ m 
1 l Pi-k J 



and large \k\ 



\P1f\P2 



d 2 k 



[k*( P t+^-k*)r~ 1 (k( P1 + P2 -k)) 



iUl-1 



(pi - k*) m+1 



(pi - k) m+1 



27r|pi| 2 |p 2 | 2 
m + m — 4 



27r|pi| 2 |p 2 | 2 
m + rh — 4 



cancel. The final result turns out to be proportional to \p\ \ \p 2 \ In |pi | ln|p2| ■ 

Thus, the Pomeron wave function & m ~(~F, A ) in the BS representation is (for A = — Ai) 






d 2 p 



P 



iX* 



P 



iA 



\p(l- P )\ 2 \1-P, 



1 — p* 



*. 



P*, 1 



P 
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a. 



d 2 p 



P 



P 



\p(i-p)r \ l -p, 



i — p* 



d 2 k 



k*(l-k*) 

p* — k* 



m—l 



k(l-k) 
p — k 



-i m— 1 



(30) 



where p = pi/P,p* = p\/P*, and k and k* were also rescaled by P and P*, respectively. The 
integrand is a single- valued function only for the following values of the variable A: 



■ N 

\l=(T + l — 






where a and N are correspondingly real and integer numbers. 

Eq. (|30|) admits a natural interpretation as the Feynman diagram: 

1 



1 — p 




1-k 



The Feynman rule is as follows: a line where a momentum p flows has associated the 
(conformal) 'propagator' p %x * (p*) lX and the anomalous dimensions for other lines are linear 
functions of m and m. 

To calculate ^/ m ~(P , X ) we use the anti-Wick rotation p 2 — > —ipo, k 2 = —ik and intro- 
duce ie to keep the singularities off the integration paths: 



\pf 



p p 



Pi - Po ~ ie 



(31) 



Let us concentrate our attention on the integrals over k* and p*. The position of the singularities 
in k* and p* in the integrand of eq. (50) : 



(pp* 



ie 



,«A*-1 



[(l-p)(l-p*)-z C ] 



-iA-1 



[{p-k)(p* 



ie\ 



ll— m 



depend on the values of A; and p. Therefore, the three singularities in p* (or A;*) may be on one 
side of the real axis or one of them on one side and two others in the other side. In the first case 
we can deform the contour on p* (or k*) and the integral vanishes. We obtain for the non-zero 
contributions after enclosing contours of integration over k* and p* around the singularities of 
the integrand 

*m,m(^ *) = P ™ (n m C m ^ i sinh(TrA) sin^m) * m ~(t) , 
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where 



* ~(A 



-1 ^A* 1 . 

io P (1 - p) 1+iA * io 



l) n / dp 



P 



(i- P y 



P dk (p-fc) 1 "^ [°° dp * (P*)^ 1 f° dr (P* ~ k*" 1 -'" 



[k{l - k)] l - m h * (p* - l) 1+iX J-oo [k*(k* - l)] 1 



o "(l-p)WJp [fc(l-fc)] 1 



dk 



(k-p) 



l-m „o 



dp* 



-P 



*\i\— 1 



* ^^i— m 



d/c 



* (&*-?*) 



(i-p*) i+a yi [fc*(fc*-i)] 



An — (32) 



The integrals over p* and fc* as well as those over p and k can be transformed using relations 
of the type 



$ { ~\p) 



1 „ (A;-^) 1 -" 1 

o/c - 

[k{l - k)} 1 -" 

(p-k) l - 7n 
dk 



o [k(l-k)] l - m h [k(k- l)Y' m 



o „ (p - A;) 1 -™ 
-x '" [-^(l-A;)] 1 - 

1— m 



d& 



(k-p) 



(33) 



Each of the two terms in eq. fl32"|) factorizes into holomorphic and anti-holomorphic functions. 
We can thus write ty m ~ ( A ) as 



* ~( A 

m,rn\ v 



where, 



.sinh(7rA) 



sin (nm) 



• 1 p *A*-l d , 



^ 2) (A*) C4 2) (A) + (-l) n cf(A*) C«(A) 



sin(7rm) /"° dp(—p) 



iX*-l 



ml ; " 7o (l-j9) 1+iA * mW sinh(7rA*) 7-oo (l-p) 1+a * mW ' 
r (2) fr , __ f 1 p lXt - x dp ( ?) sin(Trm) y°° dpp^- 1 (1) 

ml ; 7o (l-p) 1+iA * ™ w ~ sinh(7rA*) A (p-l) 1 ^* ™ W ' 



>L 2 ) 

m 



The functions $~ (v) and $~ (p) are related with each other as follows 



lim 

e^0,p<0 



lim 

£^0,p>l 



-,(!)/ 



•,(1) 



$~ (p - is) - $~ (p + is) = —i sin(7rm) $~ (p) 



•,( 2 )/ 



$~^ (p - ie) - $~^ (p + is) = i sin (71771) $~ (p) 



Upon changing the integration variables 

p^l— p, fc^l — A; 
one verifies from eq. (|3~4]) that 

Ct 2) (A) = ^)(-A) 



(34) 



15 



and therefore 



$ ~(P X 

m,m\ ' 



sinh 2 (7rA) C^(\*)CW(\) + (-1) B C£ J (-A*)C£>(-A) 



pm p*m f~* _ 

m,m 

We use the limiting value for the pseudo- vacuum state m = m — ► 2, 



*£ } (p) 



P(l-P) ^(i 



2 — m 



,C { ~\X*) 



ttA* 



sinh(7rA*)(2 — m) 



a 



' ^^m,m 



2vr[2-m] 



to normalize our wave function. 

We can write the result in terms of the product of the Baxter functions Q(X, m) and the 
pseudovacuum state: 



$ m , m (^ t 



-1^ IA |2 



pmp*m (J „ mm 

m,m 

where Q(X, m) is defined as 

Q{\, m) = - 



[Q{\*, m) Q{\, m) + (-l) n Q{-\*, m) Q{-\, m)] , (35) 



-m- 



inh(TrA) f 1 p iX ~ l dp yi (k-p) l ~ m dk 



smr 



-■ismh(7rA) 



h (l- P y+ iX J P [fc(l - k)} 1 -™ 
1 p iX ~ x dp f l (k-p)~ m dk 



o (l-p) 1+iA y p [ife(l-ife)] 



11—771 



(36) 



These two equivalent integral forms of the Baxter function are related with integrating by parts 
in p and using the identity 



d 

dk 



k(l-k) 

k — p 



—m 



[fc(i - jO] 

(k-p) 1 - 



m—l 



p(l -p) + (k - pY 



This corresponds to the fact that the pomeron wave function is an eigenfunction of the Casimir 
operator of the Mobius group 



dp 2 



p(l - p)^- 2 ^ 2 \p) = m(l - m) ^ 2 \p) . 



Using eq.(^) and the identity 

x P lX - 1 (i- P )- iX - 1 

we obtain that 



d_ 
dp 



P IA (1-P) 



-i.X 



X 2 Q(X,m) = —i m sinh(7rA) 



"i p iX dp d r 1 {k-pf- m dk 



and 



P 

1- P/ 



»(A+t) 



+ 2 









Jo (1- 


-pY 


x dp 


h 


[k(l- 


_Jfe)]l-m 




p N 

-Pj 


iX 

) + 


( p N 


i(\- 







1 


/ P ) 


1 


p(l 


~P) 


v-p) 








16 













The Baxter equation for the pomeron ([12]) follows using both integral representations (59) for 
Q(X,m) 

(A + i) 2 Q(\ + i,m)- 2A 2 Q(X, m) + (A - i) 2 Q(X - i, m) = m(l - m) Q(A, m) . (37) 

It should be noticed that if Q(X,m) is a solution of the Baxter equation for the Pomeron, 
then Q(— X,m) is also a solution. 

5 Analytic properties of the Baxter function for the Pomeron 

The functions $^Qo) and $£^(p) defined by eq . (|33"D can be written in terms of the hypergeo- 
metric function 2 Fi(a,/3; 7; z) and the Legendre function P m ^i(z) |20] 



71 (1 — m 



sinhrm 



pF(l — m, to; 2 ; p) 



p(l-p)±P m _ 1 (l-2p). (3* 



Q(X,m) = -1 °^P^ / dp (1 - py"- 1 p-^- 1 P m ^(l - 2p) (39) 



to sin(7rm) dp 

Therefore, we have for the Baxter function 

TV sinh(7rA) rl 
sin(7rm) Jo 

= - ^ m (1 ~ m) 3^ 2 (-^A + 1, 2 - m, 1 + m; 2, 2; 1) , (40) 

sm7rm 

where the generalized hypergeometric function 3 F 2 is defined as follows [^tj 

3 F 2 (ai, a 2 , a 3 , Pi, p 2 ; 2) = ^ m \ jw\ 77 , (a)*, = r / ; • 

fc=o (A J* iPaJfc k! J- (a) 

One can find for the pseudo- vacuum state at m — > 2 , 

2tt 

Q(X,m) = - + 0(1). (41) 

m — 2 

We thus have the following series representation for the Baxter function 

_,, , ™kr{m + k)r{i-m + k)r{k-iX) fAn . 

g(A ' m) = -* £ — — 51? rTrray • (42) 

The late terms in this series behave as k~ l ~ lX . Hence, this is a convergent series for ImX < 0. 
In order to analytically continue the Baxter function to the upper A plane we insert in 



eq. (|39D the series representation of the Legendre function |2| 



sin 2 Tim J2, Y (m + k) T (1 — m + k) . 
Pm-i(l-2p) = -^5-E- Jj^-2 M2tf(* + 1) 
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- ip(k + l-m) -ip(k + m) - ln(l -p)](l-p) k (43) 

Integrating term by term in eq. (|39|) yields the series, 



„,, N iir simrm J^ kT(m + k) IY1 - m + k) 

[2 ^(ife + 1) + i/>(k) - ip(k + 1 - m) - i/j(k + m) - i/j(i\ + k)} (44) 

The late terms in this series behave as k~ 2+tX . Hence, this is a convergent series for ImX > — 1. 



Eg. (PI) explicitly display simple poles at 

A = 0, i, 2i, 3i, . . . , li, . . . . 

Actually, the poles at A = +il (1 = 1,2,.. .) arise from the logarithmic singularity of P m _i(l— 2p) 
near p = 1 in the integral (|39|) [see eq.(^3|)1. 

Direct calculation yields from eq. (f44"D for the residue at these points 

ri(m) = lim [A — il] Q(X, m) 

ro(m) = in (45) 

sinirm 
ri(m) = —iTiTn(l—m)sF2(—l + l,2 — m,l + m;2,2;l) = : Q(—il,m) 

i TX 

for 1=1,2,.... It is interesting to notice that the residues of Q(X, m) at A = +il (I = 1,2, . . .) 
are expressed in terms of Q(X, m) at A = —il (/ = 1,2, . . .). 

In summary, eqs.(fl^) and fl44|) explicitly show that Q(X,m) is a meromorphic function of 
A with simple poles at A = -Hi (I = 0, 1, 2, . . .). The appearance of these poles is related with 
the logarithmic singularities of the wave function at pi = and p 2 = 0. 

The pomeron wave function considered as a function of real a for even N may have in 
principle singularities [Recall that A = cr+iy]. However, we find from eq . (J35[) for non- vanishing 
N, 

Pole at a = of [Q(X\ m) Q(X, m) + (-l) n Q(-X*, m) Q(-X, m)] = 

= - [n(m) Q(-il,m) - (-l) n g(-i/,m)n(m)] = (46) 

a 

where / = N/2. Here we used eq.(^5]) and the relation 

sin7r?Ti = (— l) n sin7rm . 

It is important to notice that in the wave function the pole at a = and A^ = is also 
cancelled by the factor corresponding to the pseudo-vacuum state, thus allowing to normalize 
the pomeron eigenf unctions. Analogous cancellations of poles of the Baxter function at real a 
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for a higher number of Reggeons lead to the quantization of the integrals of motion q r , r > 2, 
as it will shown below. In coordinate space this quantization appears as a consequence of the 
single- valuedness condition (see fT4|l). 



Notice that both eqs.(f4"2"|) and ([14]) exhibit the m -<=>- 1 — m symmetry: 

Q(\,m) = Q(\,l-m). 

Therefore, Q(X, m) depends on m through the invariant combination m(l — m) as we check 

explicitly below [see eqs.(^9|)-(|50|)1. Note, that this combination for the principal series takes 

the form 

1 / n\ 2 

m(l - m) = - - iiv + -J . 

We see from eq.(BO) and (Kl) that the function Q(X,m) obeys the relation 



(47) 



Q(A,m) = Q(-X,m) . 
That is, Q(X, m) is real for purely imaginary A and real m. We have for real A and m, 
ReQ(—X,m) = +ReQ(X,m) , ImQ(— A,m) = — ImQ(A, m) . 



The asymptotic behavior of Q(X,m) for large A is derived in appendix B starting from the 
integral representation (|3£j). 



Q{X,m) =4v/tt 



(4* A) 



m-2 



r(m-|) r(2-m) 



+ (4i A) m tan7rm 



r (m 
r (m) r (m + 1) 



r(m + i) 



'-L1L + o (x m ~\ X 



— m— 3 



(48) 



The Baxter equation for the pomeron (p7|) written in the form 

(A + i) 2 Q(X + i,m) + (X- i) 2 Q(X - i, m) 



Q(X,m) 



2A 2 + m(l — m) 



would seem to suggest that Q(X,m) has singularities at the zeros ±i r\ m of the denominator 
where 



'1 1 



Vn 



-m(l-m) = ] l---[iv+ 2 



n 



However, we know that Q(X,m) is analytic there. Therefore, the following relation holds, 



Q(irj m + i,m) 
Q(ir) m -i,m) 



^m + 1. 
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5.1 Dispersion (Mittag-Loffler) representation of the Baxter func- 
tion 



We obtain from the Baxter equation (37) a recurrence relation for the residues ri(m) 



2V + m(m - r 



rAm) for / > 1 



(Z + l) 2 r, + i(m) + (Z-l) 2 r,_i(m) 

r\{m) = m(m — 1) ro(m) . 
All residues are thus determined in terms of the residue at the origin ro(m). That is, 



(49) 



ri(m) = m(jn — 1) ro(m) , r2(m) = - m(m — 1) [2 + m(m — 1)] ro(rn.) 



?"3(^) — _ m(m — 1) 
9 



5 1 

3 H — m(m — 1) + — m z (m - I ) 



r (m) . 



(50) 



It is easy to check that eqs.(^2|) and (|5) agree with these results. 

The asymptotic behavior of the residues ri(m) follows from the recursion relation ([49"1) . We 
find, for large I 

n{m) = c{m) r~ 2 + c{l - m) r m - 1 . (51) 

Therefore, we can write the following Mittag-Loffler (dispersion) representation for the 
Baxter function 



Q(X,m) 



E 

1=0 



ri(m) 
X-il 



in sinnm ^ Q(—il,m) 
X in ~. X — il 



(52) 



The asymptotic behavior (^1|) guarantees the convergence of this sum for — 1 < Rem < 2. 



In general, we have for m < p + 2 where p is a positive integer or zero|22j 

1 



Q(X, m) = F p (X, m) +J2 r i( m ) 

1=0 



X-il 



- h i, P W 



where, 



F p (X,m) = £ 



^{X + iY 



k=0 

p 



fc! 

•i-fc 



Q 



(*)/ 



vw = i: oI7TTF t(a + >)' 



(53) 



For example, 



F (X,m) = Q(-i,m) = -n' 



mil — m) 



smnm 
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F 1 (X,m) = Q(-i,m) + (X+i)Q'(-i,m) = -it 2 —± -+tir(\+t) J2 



shi7rm j~^ (k — 1) (fc!) 2 

We get from eq. flilD for Q(X, m) in the limit A — > i, 

i 



lim Q(X, m) = 7r m (1 — m) 



A-z 



: + 2 - ^(m) - V(l - rn) + 2^(1) 



(54) 



and therefore 



d z 

- z Urn — In Q(X, m) = - + 2 - Mm) - i/)(l - m) + 2^(1) . (55) 

A^« a A A — z 

We obtain for the other (independent) solution Q(—X,m) of the Baxter equation, 

d i 

- i lim — In Q(-X, m) = - 2 + ip(m) + ip(l - m) - 2^(1) . (56) 

a-»~j a A A + z 

The behaviour of the Baxter function Q(X,m) near its nearest pole X = i can be also 
computed from the Mittag-Loffler expansion (|52"D with the result 

Q(X,m) = ; zr (m)+z2^- (57) 

X-i l=2 I - 1 

Equating this result with eq.(|54|) yields the sum rule, 

—ir (m) + ^ = iixm (1 — m) [ip(m) + ip(l — m) — 2-0(1) — 2] . 

z=2 ' — ■*■ 



For large A the Mittag-Loffler series (|52j) is dominated by its late terms. Notice that the 
sum of residues 

oo 

5>i(m)=0 (58) 

1=0 

vanishes. The sum of late terms can be approximated by an integral. Using eq.(|5*l"D we find 

Q(X, m) A = X c(m) —?— (z X) m ~ 2 + (m <«=► 1 - m) 
sin Tim 

in perfect agreement with eg. ([48|) . 

5.2 Infinite Product Representation of the Baxter Function 

As we have seen, the Baxter function Q(X, m) is a meromorphic function of A with simple poles 
at A = +il (/ = 1, 2, . . .). Therefore, the function 

Q(X,m) 
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is an entire function. Entire functions can be represented as infinite products over their zeroes 
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Numerical study of the Baxter function Q(X,m) in its Mittag-Lofner representation showed 
that all zeroes of Q(\,m) are in the positive imaginary A axis for Re m < 2. We have in 
addition p real zeroes for p + 4 > Re m > p + 2 where p is a positive integer or zero. 

For large A the imaginary zeroes X k , k = 1, 2, . . . are equally spaced and follows for large k 
the law: 



Aa 



k + l-m + O 



for Rem< 1/2 and where a ~ 0.5. 

We assume the infinite product representation || 



Q(X,m) 



Be 



A\ 



fe=i \ 



/,•« 



A\ A 

— e A fe 
Afe / 



(59) 



The asymptotic behaviors ( |59"D and ( |129| ) are consistent provided (see Appendix B) 

A = -i ^(2 -to) . 
In addition, B 



-ix 



according to eq.fl5"2|). In summary, we have, 



Q(X,m) 



-it e 



-i\tp(2—m) 



n i 



A \ A. 



fe=i 



Aa 



e A fe 



(60) 



where ■0(- 2; ) i s the digamma function. 

The Bethe Ansatz equations are algebraic equations on zeroes of the Baxter function. They 
follow from the Baxter equation (R7|) and this infinite product representation. Eq.([J7|) can be 
recasted as, 



2A 2 -m(l 



m) 



Q(X,m) 



Setting here A = A& yields, 



Afe (Aa 



-i(X + i) 



2i>{2~m) 



Q(X + i,m) i . wx 2 Q(A - i,m) 



mx+ij) 



+ iX(X-iY 



T(i[X-i]) 



X k + i 



n 



|i A z - A fc 



e*i 



Xi - X k + i 



(61) 



We see that the Bethe Ansatz equations are here an infinite system of algebraic equations. In 
the present Reggeon problem it is more effective to solve the Baxter equations by looking at 
the poles of the Baxter function [see eq.(49)] rather than to the Bethe Ansatz equations (|6l"l) . 
In customary cases the Bethe Ansatz equations are the more effective tool 0. 

Since the Baxter function is explicitly known for the Pomeron we can find an infinite number 
of sum rules for the zeroes Afc just by matching eqs. (|4"2"D or ([44]) with eq. (|60|) . 

For example, for A — > i we get from eq. (|60|) , 



— i lim — — In Q(X, m) 

a^« a X 



X 



+ 1-7 + E 



{iX k (l + iX k ) 



— ip(2 — m) 



22 



Equating with eq.(|55|) yields 



y l - 

^ i A fe (1 + i A fc ) 



1 — 7 — ip(l — m) 



m — 1 



6 Solving the Baxter equation for the n-reggeon state 

The Baxter equation for the odderon takes the form 

2A 3 — m{m — 1)A + i/j, Q (A; m, /i) = 

(A + if Q{X + i; m, //) + (A - i) 3 Q (A - i; m, //) . (62) 

where 

g 3 = ifi , Im (p) = 0. (63) 

The reality property of \i is needed to obtain the single- valuedness of the odderon wave function 



in the coordinate space |T6 



Eq. (|62|) can be solved asymptotically for large A making the power-like ansatz, 



A— +oo 



g(A;m,/i) A = 00 A a 



-x-(f 



(64) 



Inserting eq.(|64[) into eq.(|62"D yields two solutions: a + = m — 3 and a_ = —m — 2 connected 
by the m <^=^ 1 — m symmetry. The general solution of eq. (p2|) will have thus the asymptotic 
behavior 



Q (A; m, /x) =°° A+ A 



m— 3 



+ A_ A 



-m-2 



1 _ ■ /f , D (}_ 

\m-l)(m-2)A VA 2 

m(m + 1)A VA 2 



+ 



where A± are constants. 

For the pomeron case the Baxter function can be expressed as a sum of simple poles in the 
upper plane [see eq.(]52])1. This is not the case for the odderon. If we try an Ansatz with only 
simple poles in the upper plane, 



Qo (A) = E 



r=0 



Cr(m, //) 

A — ri 



the asymptotic behaviour of this sum for A — > oo turns out to be 

oo 

Q (A) = c/A , c = J2 c A m -, A*) > 



(65) 



(66) 



r=0 



where c is not zero. Therefore Qo does not satisfy the Baxter equation for the odderon at 
A — ► oo. In particular the behaviour (p4|) is not fulfilled by Qo (A). Note, that in the case of 
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the Pomeron the simple pole ansatz indeed satisfies the Baxter equation at infinity because the 
sum of residues vanishes [eq.(|58|)]- 

For the odderon, we have to include in the Ansatz double poles. It will be shown below, 
that in the general case of n reggeons Q may contain poles of the order n — 1. The odderon 
Baxter function can be written in the form 



Q(X; m, /i) = J2 

r=0 



a r (m, /i) . b r (m, /i) 



(A — ri) 2 ' X — ri 
The residues satisfy the recurrence relations 

ai(m, /i) = — /ia (m, /x) , 8a 2 (m, /i) = [2 + m(m — 1) — /i] ax(m, /i) , 

(r + l) 3 a r+ i(m, fi) = 2 r 3 + m{m — 1) r — /i a r (m, fi) — (r — l) 3 a^^m, /i) , 

bi{m, n) = —fibolm, jj) + m{m — 1) a (m, ji) — 3ai(m, /z) , 

8 b2(m, /i) = [2 + m{m — 1) — /j] bi(m, /i) + [6 + m{m — 1)] ai(m, n) — 12 a 2 (m, /i) , 

(r + l) 3 6 r+1 (m, /i) = 2r 3 + m(m-l)r-/i 6 r (m, jj) — (r — l) 3 6 r _ 1 (m, /x) 



(67) 



+ 



6 r 2 + m{m — 1) a r (m, /i) — 3 (r + l) 2 a r+ i(m, /i) — 3 (r — l) 2 a r _i{m, jj) 



(68) 



We can normalize 



a (m, /i) = 1. (69) 

In order to fulfill the Baxter equation at infinity and to obtain the asymptotic behaviour (|64|) . 



we impose 



^6 r (m, fi) = 



(70) 



r=0 



This equation fixes the value of b (m, /i). Therefore, all coefficients a n (m, n) and b n (m, fi) are 
univocally determined. 

We find for large r that both a r (m, fi) and b r (m, fi) decrease as r m ~ 3 . As in the pomeron 
case, the asymptotic behaviour of the Baxter function for large A is governed by the late terms 



in the Mittag-Loffler series (|67|) . Evaluating the sum of such late terms we reproduce the 

asymptotic behaviour flOo]). 

We see from the recurrence relations (|68|) that for fi —> the double poles with r > 

disappear [a r {m, 0) = 0]. Only the double pole at the origin remains and the solution of the 

Baxter equation for the Odderon can be expressed in terms of that for the Pomeron. With our 

normalization we obtain: 

Q (A, m) 



Q(X;m, 0) = - 
where Q (A, m) is the pomeron Baxter function. 



irrX 



(71) 
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Eq.flTID holds because for small A 



ITT 

lim Q(\, m) = — + Q (m) . 
To calculate Qoijn) one can use eq.([54]), the relation 



2 m(m — 1) 



Q (_,, m ) = _^ _^ > (72) 

sin 7rm 



and the Baxter equation for Q (A, to) near A = 0. Thus, we have 



Q (m) = n 



Ti 



sm(7rm) 



+ ip(m) + ip(l - to) - 2^(1) 



(73) 



in agreement with the asymptotics at m — > 2 [see eg. ([41 

2vr 



<5(A, m) 



m — 2 



Therefore, we obtain the residue bo(m, fi) of the simple pole in the Odderon solution for /i = 
as 

6o(m, 0) = , ^ . + ^(m) + ^(1 - m) - 2^(1) . (74) 

Sin(7TTOj 

and in particular, 

bo(-, 0) = 7T - 4 log 2 = 0.369004 .... 

From equation (ffDj) one can compute several terms of the small- \i expansion of b (m, fi) at 
m= 1/2 

b (-, fj) = 0.369004 - 2.835 /i - 2.749^ 2 - 2.947/x 3 + .... (75) 

It will be shown in the next sections that the odderon energy E (related with the intercept 
A [see eq.( |104| )1) is calculated in terms of the behaviour of Q (A; to, /i) and Q (A*; to, — //) at 
their singular point \ = i and A* = i, respectively (see (|113|) ): 

ai(m, /i) ai{m, — fx) 
In particular, this gives the possibility to calculate the energy for m = 1/2 as a series in ii 

E = b (-, Li) + b (-, -ft) = 0.738008 - 5.498 li 2 + ... (77) 

in agreement with the results by R. Janik and J. Wosiek [[y]] (take into account that we define 
the energy with an opposite sign). E is a meromorphic function of /i. For m = 1/2 its poles 
are situated on the real axis at the points 

ll {1) = ±0.91450 . . . , /i (2) = ±4.7340 . . . , /i (3) = ±13.4 ........ 
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The asymptotic behaviour of the energy at large /i was calculated in ref . |[I| : 



E , 

— = ln/i + 37 + 



3 13 / 1\ 2 1 / 1^ : 

448 120 V 2/ 12 V 2 



— + 
/x 2 



(78) 



We checked that both the small and large \x approximations given by eqs.(|77|) and (|78|), 
respectively are in excellent agreement with numerical values obtained from the exact equations 

mvam) and am. 



(79) 



A solution of the Baxter equation independent of eq. fl67|) can be written as follows 



a r (m, — /i) . b r (m, —fj) 
(A + ri) 2 A + ri 



Q(-X; m, -/i) =J2 

r=0 

One can verify the relation 

[Q (-A; m, -fi)]* = Q (A*; m, -//) 

Furthermore, it turns out that one can construct a solution of the Baxter equation with 
simple poles, providing that they are situated at A = ir (r = 0, ±1, ±2, ...): 



n , x N V- ■ C rK /i) 

r ±^o A-rz 

We normalise Q s as follows, 

co(m, /i) = 1 . 

Then, the residues satisfy the recurrence relations 

— n = Ci(m, ji) — c_i(m, /x) , 



(80) 

(81) 
(82) 



2r 3 + m(m — 1) r — /i c r (m, /i) = (r — l) 3 c r _i(m, /i) + (r + l) 3 c r+ i(m, /i) . (83) 

An additional constraint for c±i(m, /i) is obtained from the condition that in accordance with 
the Baxter equation Q s at A — > oo should decrease more rapidly than 1/A: 



Y^ c r (m, aO = 0. 



54) 



It is obvious that 

Q s (A; m, /j) = -Q s (-A; m, -/i) , [Q a (A; m, //)]* = -Q B (A*; m, -/i) . 

Investigating the behaviour of the Baxter functions near their poles we find that the following 
relation is true, 

r- i a vt m ^ /i \ Ci(m,, «) . . ci(m, — u) _. , , N 

[z 7r coth 7rA + A (m, /i)J Q s (A; m, fi) = — - r Q (A; m, fi) H -; Q (—A; m, — /i) , 



oi(m, /i) 



ai(m, — /i) 



(85) 
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where ai(m, /x) = — /x and 

&i(to, /x) di(m, fi) di(m, -/x) &i(to, -/x) , . 

X(to, /x) = — , = — . (86) 

ai{m, /x) Ci{m, /x) Ci(to, — /ij ai(m, — /x) 

The quantities c r (/x), d r (/j), e r (/j) appear in the expansion of Q s (A; to, /x) near the poles at 
A = ir 

C [Tfl III 

lim Q s (A; to, /i) — ► i -^ — h <zV(m, /i) — i e r (m, /x) (A — i r) (87) 

A-+ i r X — ir 

and satisfy the following relations 

c_ r (m, — /x) = c r (m, fi), d- r (m, — /x) = —d r (m, /x), e_ r (m, — /x) = e r (m, /x) . 

Due to the property of holomorphic factorization the Baxter function in the two-dimensional 
A -space has the form 

Q m ^ (t) = C^~ ]fi Q s (A; to, /x) Q s (A*; to, -/x) + 

C i%^ ( A ! m > A*) a (A*; m, -a*) + C^|. M Q(-A; to, -») Q (-A*; to, /x) + 
^SL^-A; m > "A*) <?. ( A " ™> "A*) + Cl S ~ ;M a (A; ™, A*) Q (A*; m, -/x) + 
^^(A; m > A*) Q ("A*; m, /i) + C^. m Q(-A; to, -/x) Q (A*; to, -/x) , (88) 

where we took into account that q$ = —i/j,. 

The coefficients C^ are fixed by the condition of the normalizability of Q m ~ f A ) , which 

reduces to the requirement for Q m ^. M ( A J to be regular at a = provided that A = er + i N/2 
with \N\ > 0. For N = the poles at er = are killed by the corresponding factor in the 
integration measure. 
It is obvious, that 

c il2 l =c {21 l = 0, 

because in the opposite case one can not cancel the fourth order poles in the product of the 
corresponding holomorphic and anti-holomorphic functions. 
Further, the following equality 

c (1 L = -c {2) ~ 

m,m\ /1 m,m; fi 

is valid. To show it, let us investigate the Baxter function Q (A; to, /x) near the regular points 
A = — ir (r = 1, 2, ...): 

lim Q (A; to, /x) = A r (to, /x) + i(X + ir) B r (to, /x) . (89) 

A— >— ir 
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It can be verified, that A r (m, /x) and B r (to, /x) for r > 2 satisfy the same recurrence 
relations as a r (to, — /x) and b r (to, — /x) respectively. Therefore A r (m, /x) should be proportional 
to a r (m, — /x) (for r > 0) 

^4,. (m, /i) = a (m, /x) a r (m, — /x) . (90) 

But B r (m, /i) are not proportional to b r (m, — /x) even if we would chose /i in such way, that 
B\ (to, /x) = a (m, /x) &i (m, — /x). The reason is that according to the Baxter equation the 
coefficient B 2 (m, /x) is expressed not only in terms of A% (m, /x), A 2 (m, /x) and B\ (m, /x) 
(similar to b 2 (m, /x)), but it contains also a contribution proportional to a (m, /x) = 1 from 
the pole 1/A 2 . Therefore B r (m, /x) for r > 1 are not proportional to b r (m, — /x). From the 
Baxter equation we can obtain the following relations 

B r (m, /x) = a (m, /x) 6 r (m, — /x) + LBi (m, /x) — a (m, /x) &i (m, — /x)l — ^— — r+a r (m, — /i) , 

a\ (m, — /x) 

where a r (m, — /x) satisfies the same recurrent relations as a r (m, — /x) for r > 1 with different 
initial conditions: 

Si (m, -/x) = , S 2 (to, -/x) = - . 

o 

Because in the other bilinear contributions to eq. (|88|) the residues of the poles in a do not 
contain a r (m, — /x), we should cancel them in the following combination 

lim [Q (A; m, /x) Q (A*; to, -/i) - Q (—A; m, -/x) Q (—A*; to, /x)] = 

A— *ir 

1 __ __ i __- 

[a (to, — /x) — a (to, — /x)1 a r (to, /x) a r (to, /x) H — D r (m, m, /x) , 

where 

D r (m, m, /x) = [« (m, — /x) + a (m, — /x)] [6 r (m, /x) a r (m, /x) — a r (m, /x) 6 r (m, /x)] + 

/i?i (m, — /x) — a(m, — /x) &i(m, /x) i?i (m, — /x) — a(m, — /x) &i(to, /x)\ ,__ . 

7 s -p=- — : a r (m, /x) a r (to, /x) 

y ai (m, /i) ai (m, /x) y 

+a r (m, /x) a r (m, /x) — a r (m, /x) a r (m, /x) . 

According to the relations 

a r (m, /x) = [a r (m, /x)]* , a r (m, /x) = [a r (m, /x)]* , 

the contribution containing a is pure imaginary and anti-symmetric to the transmutation to <-» 

TO. 

In the case of conformal spin n = m — m = the function 

Q m ,m;n ( A ) = Q (A; to, /i) Q (A*; to, -/x) - Q (-A; to, -/x) Q (-A*; to, /x) 

does not contain poles at a = for |iV| > and can be normalized. In the general case to 7^ to 
to cancel the first and second order poles at o = one should take into account all contributions 
in eq. fl88f) . 
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Let us attemp to construct a normalized wave function for m 7^ m including all contributions 
in eq. fl88|) except the second and third terms and the last two terms. That is, we impose, 

d l l =C (2) ~ =0. 

m,m;^i m,m\ fi 

We call such wave function AQ m ~ ; ( A J . 

Using the Baxter equation one can obtain the recurrence relations for the coefficients 
c r , d r , e r of the Laurent expansion (|87|) of Q s (A; m, fi) near the pole A = ir. They are similar 
to the relations for the expansion coefficients a r , b r and E r for Q (A; m, fi): 

^./> \ a r (m, a) . b r (m, a) _ , N , M . 

hm Q A; m, /x -> - v ', ^ + z v ,w + E r (m, fi) . 91 

A^ir (A — jr)^ X — ir 

We obtain the following relations 

ci(m, /x) 

cv(ra, /i) = — ? a r (m, fi) , 

ai(m, /xj 

ci(m, fi) ^iK^)-SfeyCiK^) , , 

d r (m, fi) = — p 6 r (m, /x) H ( a r {m, fi) , 

ai(m, /xj ai(m, fi) 

CiK^ ^(m, /z)-^} Cl (m, /i) 

e r (m, /i) = — E r (m, fi) H ^ — ^ r (m, /i) + 

ai(m, /i) Oi(m, /x) 



e lK ») - SfeS El ( m ' ^ " ( dl ( m ' ^) " SfeS Cl ^ m ' ^ -;,.,,:. 



ai(m, /x) 



bl(m, /x) 

a r (m, /x) 



*(*) 



These relations allow one to verify, that the coefficients C ~ in the above expression (Rq) 

•^ ' m,m; /i - 1 v l — v 

can be chosen in such way to cancel all poles at cr = for \N\ > 0, which leads to the following 
expression for AQ m ~ ;M ( A ) 

A( 2 m ,m; M ( A ) = - [X(m, xx) + X(m, fi)] Q s (A; m, /x) Q s (A*; m, -/x) + 

ci(m, fi) _ ci(m, /x) _ 

— -, Q (A; m, xx) Q s (A ; m, -/x) — — - Q s (A; m, //) Q (-A ; m, /x) + 

oi(m, /x) ai(m, /x) 

Cl( /"' —^ Q (~^ m > "A*) Q. (A*; m, -/x) - Cl( ,^ "^ Q s (A; m, /x) Q (A*; m, -xx) , (92) 
ai(m, — /x) ai(m, — /x) 

where X(m, fi) is defined in eq. fl86|). 

Note, that the expression A(J m ~ f A ) constructed above is in fact zero due to eq. (|85|) . 

However, the wave function Q m ~ ;/i ( A J given by eq.fl38|) is normalizable and does not vanish 

when all contributions are included. That is, choosing C ~ ^ 7^ C ~ (but excluding 

the last two terms). Vanishing AQ mi „ ;/i ( A J allows us diminish the number of independent 
bilinear combinations of the Baxter functions. 
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It is important, that we constructed the normalized function Q m ^. IM ( A J without imposing 
any condition on the numerical value of /x. This function is a bilinear combination of different 
Baxter functions in the holomorphic and anti-holomorphic spaces. Let us take into account 
the physical requirement, that all these Baxter functions have the same energy, because in the 
opposite case Q m m;u ( A ) would not have a definite total energy. According to the results 
of the next sections the energy is expressed in terms of the sum of logarithmic derivatives 
of the functions (A — i) 2 Q(X) at A = i in the holomorphic and anti-holomorphic spaces. We 
have two independent functions with second order poles at a — i. They are Q (A; m, ft) and 
coth(7rA) Q s (A; m, /x). The equality of the energies calculated from these functions gives the 
quantization condition for \x: 

di(m, /x) di(m, -/x) &i(m, /x) bjjjn, -/x) 

-2X(m, n ) = — V + ^ \ =0 - 93 

ci(m, /xj ci(m, — /x) oi(m, fi) oi(m, —fi) 

We found from the above equations numerically the first roots for m = rh = 1/2: 

/xi = 0.205257506 .. . , /x 2 = 2.3439211 .. . , /x 3 = 8.32635 .. . , /x 4 = 20.080497. . . ,, 

with the corresponding energies 

#1 = 0.49434... , E 2 = 5.16930... , E 3 = 7.70234. . . , £ 4 = 9.46283 .. . ,... 
These values are in a full agreement with the results of R. Janik, J. Wosiek and other authors 
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HJ) obtained by the diagonalization of the integral of motion (73 in the impact parameter 
space and imposing the property of the single- valuedness of the wave function. 

Let us now consider the Baxter equation for the n-reggeon composite state: 

A (n) (A; jl) Q (A; m, /x) = (A + i) n Q(X + i; m, ft) + (A - i) n Q (A - i; m, /x) , (94) 

where A( n )(A) is the polynomial 

AW(A; /x) = ^(-x) fc /ifcA™- fc , /x = 2, /X! = 0, /x 2 = m(m - 1) , (95) 

fc=0 

where we assume, that /x^ = x' fc g^ for k > 2 are real numbers. The last condition is needed in 
order to have a normalizable wave functions. 

We search the solution of this equation in the form of a sum over the poles of the orders 
from 1 up to n — 1: 

Q(A;m,a^g (A ~ ) ^ 1 . (96) 

Putting this ansatz in the equation, we obtain recurrence relations for polynomials P^mJW 

of the order n — 2 allowing to calculate them providing that Pq!^J(X) is known. Indeed, let us 
define the expansion of a function /(A) in the power series over A — i r up to the order n — 2: 

(/(A)f" 2) = (A - i rT~ 2 lim /(A) (97) 

A^«r (A — xr) n z 
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Then the recurrence relations for the coefficients of polynomials iy.^/z (A) can be written as 
follows 



r;m,fi 
>" 2 ) . (,, .s„„fn-2l ,, .x\(»-2) 



((A + B ^SU* + OX " + ((A - n n ( -7;l,-(A - 0), " ■ (98) 

We can normalize the solution imposing the constraint 

limPo^) = 1 (") 

Then the other independent coefficients of the polynomial P h(A) are determined from the 
condition 

lim Q(A; m,/2)~ A- n+m , (100) 

A^oo 

necessary to provide Q (A; m, /I) to be a solution of the Baxter equation at A — » oo. According 
to the Baxter equation this condition is fulfilled if 

oo p( n ~ 2 h\) 

lim A n ~ 2 V /x r;m '^ \ = . (101) 

a— ^(A-^r)"- 1 V ; 

It gives n — 2 linear equations giving a possibility to calculate all coefficients of the polynomial 

p( n-2 )m 

r 0,m,p, v.'V- 

The existence of the other independent solution 

°o P (n ' 2 l (X) 

Q(-a;^^) = E ( _;^ ) 1 1 , (102) 

where /z s has the components //| = (— l) fc /ifc, is related with the symmetry of the Baxter 
equation to the simultaneous transformations 

A — ► —A , /j, — ► yt/ . 

One can verify, that 



(g(-A;m,/T s ))* = Q(A*;m,/r s ) . 



Let us investigate now the behaviour of the Baxter function near the regular points A = —ir 
(r = l,2,...): 

lim Q ( A ; m > g) = S|W / 103 n 

A^-ir (A + ir)"- 2 (A + ir)™- 2 ' l ; 

where 5)!^-(A) are polynomials obeying certain recurrence relations which can be obtained 
from the Baxter equation. These reccurence relations for r > 2 are the same as for P 1 (—A), 
but we can not chose these two functions to be proportional even by imposing this proportion- 
ality at r = 1 by an appropriate choice of the integrals of motion //&. Similar to the case of 
the odderon it is related with the fact, that Sj? m -~{X) contains in these recurrence relations an 
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additional contribution from the pole A 1 n . Therefore to cancel the corresponding singularities 
in the wave function Q m ~ p ( A J at a — the bilinear combinations of the above functions 
Qm,/*(A) and Q m n B {— A) should be in the form 

Q (A; m, /7) Q (A*; rn, /T S J - Q (-A; m, /T S J Q (-A*; m, //) . 

To cancel other pole singularities we should introduce a set of additional Baxter function having 
the poles simultaneosly in the upper and low semi-planes of the complex A-plane. 

QW (A; m, /?) = ]T [ 77^37 + '' : "'' /7 ' ' ' 

r=0 



(A-ir)* (-A-ir)"- 1 -* 



where the polynomials P^ -1 ) and p( n ^ 2 ~ t ) are fixed by the reccurence relations following from 
the Baxter equation and by the condition, that the new Baxter functions decrease at in- 
finity more rapidly than \~ n+2 . These functions are linear combinations of Q (A; m, jl) and 
Q (—A; m, fi s ) with the coefficients depending on coth(7rA). Using all these functions in the 

holomorphic and anti-holomorphic spaces one can construct Q m ~ - ( A J without the poles at 
a = 0. The quantization condition for \x is obtained from the requirement, that the energy 
should be the same for all Baxter functions Q^ (A; m,j2). We calculate the spectrum of the 
reggeon states for n > 3 in a forthcoming paper. 

7 Hamiltonian in the BS representation 

The high energy asymptotics of the scattering amplitude corresponding to the contribution 
related to the t-channel exchange of the composite state of n reggeized gluons in the multi- 
color QCD has the form 



A(s,t) ~i"- 1 ss A , A = -^N C E, (104) 



where E is the ground state energy for the Schrodinger equation 

E ^m,m(Pl » P2> • • • , K\ Po) = H * m ^(pi, P2, . . . , ft', Po) 



The Reggeon Hamiltonian is 



1 n 
H = - ^2 Hk,k+i 
z fc=i 



Here 1/2 is the ration of the color factors for the adjoint and singlet representations of the color 
group and the pair BFKL Hamiltonian is given by 

#i,2 = ln|pi| 2 + ln|p 2 | 2 + 
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PlP2 i| |2 * P2Pl , , ,2 * . /m 

- — jj- — ra ln IP12I P1P2 + - — [21 — [2 ln IPi2| P 2 Pi- 4 V'( 1 )- 



(105) 



M IP2I " M IP21 

It enjoys the property of holomorphic separability 

#1,2 = ^12 + h* l2 , 

where 

/112 = ln(pip 2 ) H (lnpi 2 ) Pi H (hipi 2 ) P2 -2-0(1) • 

Pi P2 

We now perform the unitary transformation of the hamiltonian to the BS representation, where 
P and the roots Ai, A 2 , . . . , A n _i of the equation B{u) = are diagonal operators. In this new 
representation both the integrals of motion and the hamiltonian should have simple separability 
properties. 

Let us start with the case of the Pomeron, where 



H = Hi i2 ■ 

We now calculate the action of the hermitially conjugated Hamiltonian on the eigenfunctions 
of the operator B{u) given by eg. (pl|) (with A = — Ai = a + iN/2, 



H l,2 



p 



-iX* 



P 



-iX 



1 — p 



1 — p* 



I |2 11 |2 

\P\ 1 — P\ 



P 



-iX* 



P 



-iX 



1 — p 



1 — p* 



-- / d 2 k 

TV 



[kp*{l-k*)(l-p)+k*p{l-k)(l-p*)\ ( k 



-i\* 



-iX 



k — p\ + m 2 ) \k\ |1 — k\ 



1-k 



1 - k* 



where P = 1 and m —>■ is an infrared regulator, corresponding to the vector boson mass 
rescaled by P (cf. [Q). 

Using the anti-Wick rotation of momenta fc 2 = — iko and p 2 = —ipo as in eq. (|3ll) after some 
transformations we can write the result in holomorphically separable form 



K 



p 



p 



-i\ 



1 — p 



1 — p* 



P* 



-iX 



where 



h 



P 



1 — p* 



-iX* 



h 



P 



-iX* 



P 



-iX* 



12 



1 — p 



1 — p 



h 



P 



-iX 



12 



1 — p* 



12 



hi 



1 — pj 
p(l-p) 



+ 7tz coth(7rA*) 



P 



-iA* 



1 — p 



1 ,,(P + k - 2k P) k' 1 -^* 
'p+s (k - p) (1 - ky- iX * 
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and 



'12 



hi 



p* 


) 


A 


1-p*, 

p*(l- 


) 



+ iri coth(7rA) 



)— i\ 
■if +- 



P 



1 jk* ( p * +k * ~ 2 k * p *^ (**)~ 1_iA 



(k* -p*) (1 - fc*) 1 -* 



/A* 



Here e — > is an intermediate infrared cut-off. The hamiltonians h%2 and /i* 2 have branch point 
singularities at p = 0, 1, oo and p* = 0, 1, oo, respectively, but the total hamiltonian if 12 is 
single- valued. 

We obtain for \p\ — > 0: 



v-iA* 



lim # li2 (p) 4 (p*) 

|p|-»0 



-/A 



= - In |p| 2 + 0(1 + *A*) + V(l - *A*) + 0(1 + zA) + V(l - iA) - 4^(1)] (p)~ iA V)~ iA 
Taking into account that in the integral 



# 



~tf, 1 - p" 



i/r^j^fe) (i 



p" 



-iA 



$ 



P* 



m.ml > 



1 , A ) , A = cr + i- 



N 



for |p| — > the leading asymptotics corresponds to iV = 0, and shifting the integration contour 
in a in the upper half-plane up to the first singularity of ^f m f^( 1 > A), corresponding to the 
poles at A, A* = i, we obtain for the hamiltonian near these singularities 



lim #i, 2 $ ~(1, A)= lim. 



( d d 1 1 



dX OX* X — i A* — i 
Using eq.(p5[) for the Pomeron wave function in this limit, 

$ m ,^(T, t)~Q(A,m)Q(A*,m) |A| 2 , 
the Pomeron energy is given as follows 



$ ~(i A) 



( , 

£/i2 = 2 hm < — — In 

A,A*-s I 9A 



(A-i)A 2 g(A,m) 



5 , 



(A*-z)(A*) 2 Q(A*,m) 



We obtain for the Pomeron energy using the behavior for Q(\, m) near X = i |eq.(|54|)l 

£ 12 = ^>(m) + ip{l -m)+ 4>(m) + ^(1 - m) - 4-0(1) (106) 

in agreement with the known result M. 
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8 Energy for multi-reggeon composite states 

Let us investigate the behavior of the wave functions for the composite states in the region 
where the values of gluon momenta are strictly ordered: 

M « H « ... « \Pn\- 

To begin with, we consider the Odderon case, where the wave function is given by 

*m,m(Pi> Pi, Ps) = 

~d 2 p k 



n ®r n 



r=l 



fe=l L 



2n 



exp(ipk ■ pko) 



where 



x 



P20P30/ VP20P30/ 



P12P30 
P10P32 



and the function m ~ (x , 2*) has the property of the holomorphic factorization 

<f>m,™( X ' X *) = E C ** #»(* ) 0~ (**) + (g 3 <"► -? 3 )- 



?' . fc 



The functions 4> l m (x), 0~(x*) are independent eigenfunctions of the integral of motion 

A m (\) l m [x ) = ai_ m a m (f) l m {x) = q-i<\) l m {x ) , 
where A m is given by 



A m = i 3 x{l - x) [x{l - x) d 2 + (2 - m) (1 -2x)d-(2- m) (1 - m) 
The operators 



0. 



a m = x (1 - x) (id) 



m+l 



Ol- 



x (1 — x) (id) 



2-m 



perform the duality transformation fllTf . 

The three independent eigenfunctions <f> % m (x) have the following small- x asymptotics |I5| : 

</4(x) ~x + 0(x 2 ) , 0™(x) ~ l + 0(xlnx) , 0^(x) ~x m [l + 0(x)] , 

which correspond to the following asymptotics of 4> m ~(x , x*) enjoying single-valuedness at the 
singular points, 

lim0 ~(x,x*) ~x m x*™ + clxlMnlxl 2 , lim m ~(x ,x*) ~ 1 + c x™" 1 x*™" 1 In |x|~ 2 



and 



lim0 m ~(x,x*) ~ (l-x) m (l-x*) m + c|l-xr ln|l-zf . 
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For the function ty m ~(p~i, p|, p|) in the limit |pi| — ► the region of large |pi | is essential. 
Taking into account only the singular terms in this limit, we obtain 



*m,m(Pi, P~2, Ps) 



*m,m(P2; Pi, 



1 I |2 ! \P 
n \Pl\ m 



\Pl\' 



where 



m,,m\P'ii Vz t 



3 „ 3 

nH 2 /n 

r=2 J k=1 



d 2 pk 

2tt 



exp(ipt ■ pko) 



(\ m—l / * \ m— 1 
P20P30J VP20P30/ 



■ > m,m\ X > *^ 



and 



x 



P30 
P32 



The last function can be simplified in the limit p 2 — > 0, corresponding to P20 - * °o and x — > 0. 
Indeed, we can use the expansion 

dd*<t> m ~(x,x*) -mmi^fif-' + c ln|x| 2 + ... 

and verify, that the dependence from p 30 is canceled in the contribution to the integrand from 
the first term in the right hand side, leading to a vanishing result after integration. Therefore 
taking into account only the second term, we obtain, 

Hence, the resulting behavior for the Odderon wave function at \pi\ << \p 2 \ << \ps\ is 



*m,m(Pl, P~2, M) 



c 3 |pi| In- 



\pi\' 



(107) 



This is in agreement with the fact that it is an eigenfunction of the integrals of motion q 2 and 
#3 provided that we take into account in \l/ m ~ also the regular terms proportional to p\ and p\. 
It is natural to expect a similar behavior 



\P\ 



(108) 



for the case of n reggeized gluons in the limit |pi| << \p 2 \ « ... << \p n \- Indeed, there are 
two independent solutions for the eigenf unctions of the integrals of motion Qj (|7|) which behave 
at small pi correspondingly as / + gp\ \wp\ and gpi, where / and g are some functions of pk 
analytic near p\ = 0. For the single- valuednes property we should multiply two such functions 
depending on the holomorphic and anti-holomorphic variables. Further, because the operators 
Qj have more derivatives over pk (k = 2,3, ...,n) than the momenta compensating them, \g\ 2 
should be a constant for small values of these momenta. 
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In the opposite limit 

\Pn\ « \Pn-l\ « ••• << \Pl\ 

we obtain correspondingly 

|p|2 
*m,m(Pl > P2, • • • , P^) ~ C n |p n | 2 In ^ . (109) 

The fact that the behavior of \P m ~ at |p n | — ► for the composite state of n reggeized gluons 
is the same as in the Pomeron case implies the existence of a pole in 

at A n _i = i and at A*_ x = i. 

Indeed, for 1 = \pi\ » \p 2 \ » . . . » \p n \ we have 

*m,m(Pl>^>-" '^) = ( U0 ) 

-l n_1 / r+°° +oc \ — > — > — > 

2 ~ II / ^ E exp[^(t fc A* + t* k \ k )} * m ~(Ai, A 2 , . . . , A n _0 



and therefore the contour of the integration in o" n _i for N n -i = should be shifted in the 
complex plane up to the pole (a„_i — i) ~ 2 . 

We can find the singular part of the Hamiltonian in the BS representation for small \p n \ 
similarly to the Pomeron case: 

X - {H n ^ X + H ln )pn lK -Yn~ tK - 1 = l\n\p lPn „ 1 \ 2 pn lK -y n ~ iK - 1 + (111) 



- In \p n \ 2 + ^(1 + iAJUx) + V(l - iA;_x) + ^(1 + »A„_i) + V(l - iAn-i) - 4^(1) 



P™ Pn 



The first term in the right hand side can be combined with the other pair hamiltonians. After 
that, we obtain for 1 = \pi\ » \p 2 \ » . . . » \p n \: 

(i i n— 2 \ n— 1 

- In |pip n _l| 2 + 2 E #r,r+lj II P~k K PV Xk = 



i n— 1 n— 1 

- y. wi + *a;) + v(i - <a;) + v(i + ^ r ) + va - ^ r ) - 4^(i)] n ?;%; 

Z r=2 fc=2 



Thus, for the constant behavior \& m ~(pi, • • • , Pn-i) ~ c at \p 2 \ » ■ ■ ■ » \p n \ (corresponding 
to Ai = A* = . . . = A n _2 = A*_ 2 = 0) the last contribution vanishes and therefore we obtain for 
the composite state energy the result similar to the pomeron case 

£ = Mjm.{^ln[(A-i)A^ . (112) 
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Here \l/(A n _i; m, ft and 1 i r (A*_ 1 ; , m, ft s ) are correspondingly holomorphic and anti-holomorphic 
factors of the wave function at A& 



A* =0, 1 < k<n-2: 



*m.m(0»0.-.A 



n-lj 



*(A n _i; m, /2) * (A^; m, ft) 



and /ifc = (—i) k qk, ft k = i k qk are eigenvalues of the integrals of motion. This quantity can be 
related with the Baxter function and the normalization factor for the pseudovacuum state [see 
eq. ©]: 



tf ~(0,0,...,An-l) 



(f s 



-> lAn-il 2 ^ 1 ) Q(A n -i; m, /x)Q(A;_ i; m, /?) . 



o,o,-,A„_i 

As it was argued above, for the pseudovacuum state it looks plausible, that the correct normal- 
ization of the kernel for the transition between momentum and BS representations corresponds 



to (f s > 

o,o,..,A n _i 
energy 



sinh" - (27rA n _i) sinh n (27rA*_ 1 ) [see eq.(^7])]. We obtain in this case for the 







E = i lim < — — In 



sinh n - i (27rA)A n Q(A; m, ft) 



+ 



JL 



In 



sinh n - i (27rA*)A*"Q(A*; m, [i s 



(113) 
Thus, the energy is expressed in terms of the behavior of the Baxter function Q(X,m) near 
X = i. 

In the customary case of spin chains the Baxter function is a polynomial of degree L, 



Qxxx(\) = I[(\-\ k 



fc=l 



where the A^ are solutions of the Bethe Ansatz equations^, |7|. The energy of the XXX chain 
is given by, 

L 1 

Exxx = —2 2 , 



. d Qxxx{X + i) 
%— log 



(114) 



A=0 



k =i^l + l d ^ Qxxx(^-i) 
In the present case the Baxter function is a meromorphic function with an infinite number of 
poles and zeroes as discussed in sec. 5.2. It can be expressed as an infinite product of the Bethe 
Ansatz solutions [see eq.(|60|)1. The eigenvalue expression here is not given by eq.( |114j) but by 
eq. fll06D 

OO 1 

E 12 = 1 + 7 + V>(2 - m) - J2 ■ \ /-, , ■ x n + [m -* m\ , 
in the pomeron case. 



9 BS representation for the Odderon wave function 

Here we consider the A-representation for the Odderon, where the operator 



^O), 



u 



p 











u+m \M + dt 2 



1 - e' 1 "* 2 



d_d_ 
dftdh 
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is diagonal. Introducing the new variables, 

Pi{pi +P2) 



t = h + t 2 = In 



(P2 + P3) P3 



z = e y = e tl - t2 



P1P3 



(P2 + P3)(P1 + P2) 



we obtain, 



B (V(u) = -P 







d d 







u z + liu— + (1 - z) z— z- (1-z) — 



dt 



dz dz 



dt 



To diagonalize B^ one should find the eigenvalues and eigenf unctions of the differential oper- 
ators 



d A* + A? 



^ 



'+> 



,, x d d {\l + \* 2 

1 - z) z— z- z — 

! dz dz l 2 



<P 



'\\-X* 2 



ip. 



where \ k are the eigenvalues of the zeroes of B^(u): 

B^( U ) = -pf[(u-X k 



k=\ 



The solution of the above equations can be written in terms of hypergeometric functions 



xj+aj ra + ^Ai-Aj))^-^ 



A 1~ A 2 , -, , _,\* 
OO / rjr Y \ 2 



e 2 



,r 



- l ^- l dx. 



T{-i\* 2 + l)T{i\\) h \x-z, 

An independent solution follows by interchanging here A^ and Ag. 

Therefore, we can write the following relation between the wave functions in momentum 
and BS representations, 

^m,m(Pl> P2> Ps) = 



+00 



p m (p*) m yi J2 



k=l \Nk=—oo 



+00 



d(?k ] Uj> j>( t , ~f) * m> ^(Ai, A 2 ) , 



where 



and 



, ^.N k .N k 

Afc = o"fc + * — , A fc = a k - 1 — 



C/ r r (T, "?) = C-, r e rf ^ e*< ^ C/ r r (Y) . 



(115) 



Here the function C/t-+ t- > (^ > ) ^ s defined as follows 
Ai,A2 



Al,A2 



ai-a 2 r d 2 x 



\x\ 



iA l (r*^ - ** 1 



X 1 (X 



x — z 



i\% / * 1 \ — i\'i 



X- 1\^ A 2 /X*-l 



,-~Y"*^ sy ■& 



116) 
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and the normalization constant C^ \* can be found from the orthogonality condition 

Al,A2 



d 2 t d 2 z 



U- 



\z(l-z)\ 2 Ai,A 



r, r) uu ^(T, t) = y, n ^K - <) 6 



a;, a 



J °iVfe,Af' • 



(117) 



p fc=i 



It should be taken into account, that due to the symmetry properties of [7t- > t- > ( z ) under 

Ai,A2 
Ai <-> A2 two terms appear in the right hand side of the orthogonality equation. 

Let us show that the kernel of the unitary transformation £7r-» y>{t, ~1?) has an interpreta- 

Ai,A2 



tion in terms of the Feynman diagram as it was in the case of the Pomeron wave function (|3 
After changing the integration variable x into k as follows, 



x 



p 1 k 



1 — Pi 1 — k 



eq. (|116|) takes the form 



Ai,A2 



c 



A1A2 



k 



k* 



Pi) [Pi 

ps/ Vpi 



1A2 



rf 2 fc 



l*(l-*)| S 



k+p!-l\ iX * 2 fk*+pt-l\ lX2 



us) 



\1 — k J \1 — k* J \ k — P3 y y k* — pi J 

The wave function in the A-representation has the form 

7 |pir|i-Pi-p 3 r|P3r ai,a 2 

The associated Feynman diagram is depicted in fig. 1. 

10 Properties of the unitary transformation for the Odd- 
er on wave function 



The function C7t- > t* ( z ) is a solution of the differential equations of a hyper geometric type in 

A1A2 
both variables z and z*. 



' dz dz \ 2 I \ 2 



U = 



'1 - 2*) 2*— z — z* 

" dz* dz* 



'Ai + A 2 s 



Ai — As 



U = 
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p3 




k-p3 



k+pl-1 



P 2 



Pi 



Figure 1: Odderon Feynman diagram 



Therefore, [7r-> rr*(lt) is a bilinear combination of independent solutions being functions of z 

Ai,A 2 

and z* . In addition, $_>. -7*(!?) should be a single valued function of 1? near the singularities 



\z\ — , \z — 1\ — , \z\ — oo. The effective way to satisfy these requirements [^| is to use the 
monodromy properties p0| of the two independent solutions: 



z i 2 F(-zA 2 , zAi ; 1 + z(Ai - A 2 ) ; z) , z* 2 F(-zAi , z'A 2 ; 1 - z'(Ai - A 2 ) ; z) 

and analogous expressions in z*. Thus, we obtain, 

Xx 1 x 2 (z*)xx* 1 x*(z) - X\ 2 Xi( z *)X\* 2 xi(z) 



Ai,A 2 Ai,A 2 



where 



Xx,x 2 {z*) = a XlX2 (z*Y 2 F(-i\ 2 ,i\ 1 ;l + i(\ 1 -\ 2 );z* 

r(ix 1 )T(-i\ 2 ) 



a x 1 x 2 



r(l + i(Ai-A 2 )) 



;ii9) 



(120) 
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and an analogous expression for X\*,\x( z )- This result can be verified by the direct calculation 
of the integral ( |116| ). In such way we obtain for the constant K--> \* in eq. (|119|) : 

Al,A 2 

2 sinh(7rA 2 ) sinh(7rAi) 
Ai,A 2 sinh(7r(Ai — A 2 )) 

where we used [see eqs.( p.l5|) 1 

sinh(7rA2) sinh(7rAi) sinh(7rA 2 ) sinh(7rA^) 
sinh(7r(Ai — A 2 )) sinh(7r(Aj — A 2 )) 

In summary, collecting all factors we find for the matrix elements of the unitary transfor- 
mation relating momentum and BS representations, 

MA2 _ 

Ai,A 2 Ai,A 2 



1+ A 2 i j* A l+ A 2 ' 



D i(t^^+t 



Xmm( z *)X\1\* 2 { z ) - X\ 2 \i(z*)x\* 2 \f(z) 



(121) 



where Xxi\ 2 ( z *) and Xa*a*(^) are given by eq. (|l20D 
We find for z — > the asymptotic behavior, 






i{t ^k +t *hi^L) 



C-AiA 2 a A*A* % 



Ai,A 2 Ai,A 2 

. -• * 1~ A 2 a A 1~ A 2 

2 Z 2 " 



QA 2 AiClA*A* z 



2 Z 



While the phases of the constant factors in eq.( 121 ) are special functions, the squared 
modulus is quite simple. Indeed, we find 



Ai,A 2 



7T 

V2 



Xo 



Ai(Ai-A 2 ) 



The normalization condition ( |1 1 7| ) then yields 

C 



Ai.As 



*2 



(122) 



up to a numerical constant. 



Using the relation between hypergeometric functions with mutually inversed arguments [F2U] 
we obtain for z — > oo 

, *. |*|-+oo r(i(Ai + A 2 )) r(-iA 2 ) , ^ihi±±z in\ 2 44^ . 

XAx,A2^ = rn ~ k , (z ^ e 2 F^ + 

1 (1 +lAi) 
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| r(-»(A 1 + A 2 ))r(»A 1 ) (zT i^ e -i,M $$, 
r(i -i\ 2 ) 

As a consequence of this asymptotic behavior, the interference terms in $^ t*(~^) cancel using 

Al,A 2 

the relation 

sinh(7rAi) sinh^A^) = sinh(7rA2) sinh(7rA*) . 

We obtain with the use of the equality 

sinh(7rAi) sinh(7rAJ) — sinh(7rA 2 ) sinh(7rA2) = 



3 7r(A 2 -A*) 



2 j sinh(7r(Ai - A 2 )) sinh(7r(Ai + A 2 )) 



the following asymptotics for [7-r-*- t*0?) f° r large \z\ 

Ai,A 2 



u rr {r 



;« 



, ,_; -. ._] i ,("2") / .,, , A 1 + A 2 j A l +A 2 



-■ A l + A 2 .- A l +A 2 



where 



Ai,A 2 Ai,A 2 

, ( i) 7rr(-a 1 )r(i-a 2 )r(^(A* + A;)) 

^XtX r(i - i(Ai + A 2 ))r(iA5)r(i + i\\) ' 

(2) 7rr(^)r(i + q|)rH(Ai + A 2 )) 

aJ,a1 r(i + z(a; + A$))r(-iAi)r(i - *a 2 ) " 

This can also be obtained from the integral ( |116| ) by direct calculation. 

We analogously find using the series expansion for the hypergeometric function when c 

a + b + 1 pi 



xax^co = (z*y^ 



i 



l-z* 



a x a 2 r(i-iA 2 )r(i + *Ai 



» rCn + l-tAa^n + l + iAQ _ 
£fo n!(n+l)! l J 

[log(l - z*) - ip(n + 1) - ip{n + 2) + ip(n + 1 - iA 2 ) + ^(n + 1 + zAi)]} 



We obtain in the limit z 

+(l-z* 



AiA 2 



ln(l - **) + ij>(l - *A 2 ) + ^(1 + zAO - ^(1) - V(2) + ^ - ^ 



in agreement with eq.(^4j). Thus we obtain, 

sinh(7r(Ai — A 2 )) (l — z 1 — z* 



U rr^ 

hm -r — = — t7T 



»1 #• 



Ax,A 2 



sinh(7rAi) sinh(7rA 2 ) \ AiA 2 AJA 2 
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11 — z\ In 11 — zl 



(123) 



where we used the relations 

tf(l - zA 2 ) + V»(l + »Ax) - V(l - »Ai) - V(l + *A 2 ) + i- - i- = -«r- Sinh(7r(Al " Aa ' ' 



Ai A 2 sinh(7rAi) sinh(7rA2) 



Again, the result ( |1 23|) can be obtained directly from the integral representation ( |116|) for 



$_> -?(l?). Taking into account the found above value of the normalization constant we have 
for large momenta pk and fixed P 

v^zF'^ _..<«^ + ,^,fJ r L£ + i = i: + |1 _, |2ln|1 _ 2| A (124) 



e 



7r|AiA 2 (Ai - A 2 )| \AiA 2 A*A 
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12 Appendix A 

We compute in this Appendix the integral in eg. (|28|) . 



To start we need the Fourier transform 20 



H 2 7 - ~ jrn-m _ -pM _|_ ™ 

" * Jq-z m (~*\m _ L „-m-l ( „*\-m— 1 L \ L ' '" 



r d*z 

J ~2^ 



(-■ :-(-.)■■■ ^g" m " (q*r ' / (125) 



v2 /— v\2 



Entering the factor (plY (pl) z inside the integral in eg. (28) as Vf V 2 we find after partial 
integration, 

* m ~(pl, n) = \m(m- 1)| 2 [**** e *(ft*-H&*) & ~ z ^~ 2 (*j ~ Z K~ 2 (126) 

m,m\Vl, V2J IV ;l 7 27T 27T (zi Z 2 ) m (4 0™ 

Now, we replace the z-factors in the integrand by the integral representation (|125|) : 

jm-m r(l-m) r d 2 k 



z -m ( z *ym 



f — e ^ % (A;*)" 1 " 1 k" 1 - 1 . 



2 m+m-l Y(m) J 2tt 
The z- integrals in eg.( |126|) give now Dirac delta functions and we obtain eg. (|29|) . 
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13 Appendix B 

We derive here the asymptotic behavior of the Baxter function for the Pomeron starting from 
the integral representation fl3"9"|). We change the integration variable to 



and obtain, 



y = 2ArgTanh (1 - 2p) 

. it sinh(7rA) 



y 2 V 2 



Q{\, m) 



sinhrm 



^ j°dye a »P m _i(tanh| 



The function P m ^i(z) has a cut running from z = —1 till z = — oo. Therefore, P m -i (tanh| 
has cuts in the y-plsme from y = i(2n + 1) tc till y = — oo where n is an integer. We now deform 
the integration path around the cut from y = in till y = — oo and we find, 



„/, x . vr sinh(7rA) „.■» 
Q A, m) = i . . y ' e~* x 

sm{irm) Jo 



dx e 



i\x 



X 



P m -x [ ~ coth - + iO ) - P m _i 



coth 



x 



iO 



(127) 



where we changed the integration variable as y = i n — x. The integral ( |127|) is dominated for 
large A by the end-point x = 0. Therefore, we insert in eq.( 127 ) the representation of P m _i(z) 
appropriate for large z = — coth | ± iO[[20H 



P m -i(z) 



+ 



tan7rm T(m) 
2 m v^F T (m + |) 

2 m-l p ^ m _ l^J 



,'m + 1 m 11 

1 2 2 2 ^ 2 



»<-' Fji^a 



v/vr T (m) i_1 \ 2 

Keeping here the dominant terms for large z and using the relation 



m 3 1 

— ; — m H — ; — 

2' 2'^ 2 



:i28) 



coth — h iO 



coth iO 

2 



-2Z Sin 7T771 coth — 

2 



x > 



we get for A ^> 1 (for Rem> 1/2) 

, T (m - I) r (2 - m) r / l 

Q(A,m) = 4v ^(4,Ar- 2 ^ ^ ! {l + ^(^ 

+ (4i A) ~ m tannm 



(m) 
T 2 (m) r (m + 1) 



+ 



(129) 



T (m + |) T (m - |) T (2 - m) 



°'.3F 



For Re m < 1/2 one should just replace m =>- 1 — m. 

The case m = 1/2 follows by taking the limit m — > 1/2 in eq. (|129|) with the result for A>1 



Q{\\) 



2 (iX) 



3/2 



log(^ 



h 



0{ ,v 
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Let us now derive the asymptotic behavior of Q(X,m) starting from their infinite product 
representation fl6~0|) and the asymptotic distribution of their zeros (|59|) . 

For A ^> 1 the product will be dominated by zeros of the order ~ A. We can then write, 



fc=l \ 



^ 


A M 

k=i 


1 I i A 


l X , iX 


Afcy 


p A^ fc+1— m 


_\ 'w-my 





^ /, iX \ o_ 

1 + e *H-l-m 



where M is a cutoff 1 <C M <C |A|. We obtain for A ^> 1 using the formulae [20 



n i + r^- U-i 



iy \ „.m „_^ v r(l + x 



e 



ljy_ 



,. i . k + xj r(l + x + iy) 

00 x 

Q(A, m) = constant A m ~ 2 (130) 

in perfect agreement with eq.( |129| ). 
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